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We report the remarkably robust ferroelectric state in the multiferroic compound Mni-^Zn^WCU. 
The substitution of the magnetic Mn 2+ with nonmagnetic Zn 2+ reduces the magnetic exchange and 
provides control of the various magnetic and multiferroic states of MnWC>4. Only 5 % of Zn substi- 
tution results in a complete suppression of the frustrated collinear (paraelectric) low temperature 
phase. The helical magnetic and ferroelectric phase develops as the ground state. The multiferroic 
state is stable up to a high level of substitution of more than 50 %. The magnetic, thermodynamic, 
and dielectric properties as well as the ferroelectric polarization of single crystals of Mni_ a; Zn a: W04 
are studied for different substitutions up to x=0.5. The magnetic phases have been identified in 
single crystal neutron scattering experiments. The ferroelectric polarization scales with the neutron 
intensity of the incommensurate peak of the helical phase. 

PACS numbers: 75.30.-m,75.30.Kz,75.50.Ee,77.80.-e,77.84.Bw 



I. INTRODUCTION 

Multiferroic materials have received increasing atten- 
tion in recent years because of the coexistence and mutual 
interaction of magnetic and ferroelectric orders and their 
coupling to external magnetic and electric fields.— ~— Im- 
proper ferroelectricity can be induced by certain kinds of 
inversion symmetry breaking magnetic order if the mag- 
netic moments are strongly coupled to the lattice causing 
the ionic displacements with a macroscopic electrical po- 
larization. The transverse helical magnetic structure was 
shown to be compatible by symmetry with a ferroelec- 
tric polarization arising from a third order coupling term 
of the ferroelectric and magnetic order parameters in 
the Ginzburg-Landau thermodynamic potential. 4,5 He- 
lical (non-collinear) magnetic structures have in fact 
been found in various multiferroics, for example, in 
TbMn0 3 ^ Ni 3 V 2 8 ^ and MnW0 4 ^ Different micro- 
scopic models have been recently proposed to describe 
the phase complexity and ferroelectricity in multiferroic 
manganites 3~— 

The non-collinear magnetic order observed in these 
compounds is a consequence of magnetic frustration due 
to geometric constraints or competing exchange interac- 
tions resulting in a close competition of different mag- 
netic structures that are nearly equal in energy. Most 
multiferroic materials are therefore extremely sensitive to 
small perturbations. Magnetic fields may stabilize or de- 
stroy the multiferroic state or result in the rotation of the 
ferroelectric polarization by 90° or even by 180° Zd2r— 
Similarly, the application of external pressure^— as well 
as chemical substitutions 2 ^— are a viable tool to change 
and tune the multiferroic properties of different com- 
pounds. 

MnWC>4 is also known as the mineral Hiibnerite. The 
crystallographic structure of the compound is mono- 



clinic (space group: P 2/c). Three magnetic phase 
transitions at Tat=13.5 K, T c =12.6 K, and T L =7.8 K 
separate different collinear and non-collinear magnetic 
phases^ The AF3 phase (T c <T<Tat) shows a sinu- 
soidal spin order given by the incommensurate (IC) vec- 
tor if 3 = (-0.214,1/2,0.457). The collinear Mn spins 
are confined to the a-c plane forming an angle of about 
34° with the a-axis. Below Tc the spins tilt out of 
plane and form a helical, non-collinear structure break- 
ing the spatial inversion symmetry (AF2 phase). The 
modulation vector of the spin order does not change 
at Tc, ~? > 2 = ~^3- The helical phase becomes unsta- 
ble at Tr, and the low temperature AF1 phase (ground 
state) is characterized by a frustrated tT44- s P m struc- 
ture with a the commensurate (CM) modulation vec- 
tor q 1 = ( — 1/4,1/2,1/2). Among the three magnetic 
phases only the helical AF2 phase is ferroelectric. The 
details of the magnetic structures have been revealed in 
neutron scattering experiments 1 24 i 25 The complex phase 
diagram of MnWC>4 and the multiferroic properties have 
been described extensively through magnetic, dielectric, 
and thermodynamic measurements^^ - — Novel physi- 
cal phenomena such as the coupling of ferroelectric with 
magnetic domains 3 - - or the control of the magnetic do- 
main chirality by electric fields 3 - 1 - have been reported. 

The tuning of the magnetic exchange interactions 
and, consequently, the multiferroic properties by partial 
substitution of manganese by other transition metals, 
for example iron or cobalt, was achieved recently^r— 
The magnetic phases and the multiferroic proper- 
ties are extremely sensitive to small amounts of Fc 
substitutioni 20 ' 21 ' 37 Some of the observed effects of tran- 
sition metal substitution can be qualitatively described 
by simple Hciscnbcrg models with competing interactions 
and uniaxial anisotropy ; 20 ' 38 however, a more quantita- 
tive and microscopic understanding of the complex mag- 
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FIG. 1. (Color online) Lattice constants a, b, and c of 
Mni-^Zn^WCU as function of x. 



netic interactions and the effects of Fe or Co substitu- 
tion seems to be a challenge because of the fact that Fe 
as well as Co carry their own magnetic moments and 
their exchange interactions with the Mn spins as well 
as their anisotropy parameters are not known. To over- 
come some of these problems and to simplify the physics 
of magnetic exchange in substituted MnWC>4 we decided 
to replace the Mn 2+ by non-magnetic d-metals such as 
Zn 2+ . We have therefore synthesized single crystals of 
Mni-zZn-rWC^ with x between and 0.5 and investi- 
gated their magnetic, thermodynamic, and multiferroic 
properties. 
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FIG. 2. (Color online) Inverse magnetic susceptibility of 
Mni-^Zn^WCU. The inset shows the x-dependence of the es- 
timated effective magnetic moment [i e f f and the Curie- Weiss 
temperature 6w- The dashed line is the expected y, e ff for an 
Mn spin of S=5/2. The black (solid) line is a linear fit to the 
data of <dw- 



three axis spectrometers at the High Flux Isotope Reac- 
tor at the Oak Ridge National Laboratory. The crystals 
were aligned in the scattering plane defined by the two 
orthogonal wave vectors (1,0,-2) and (0,1,0), in which the 
magnetic Bragg peaks and other structural peaks can be 
surveyed. The incident neutron energy was fixed at 14.7 
meV using pyrolytic graphite crystals as monochromator, 
analyzer, and filter. 



II. EXPERIMENTAL 

The crystals were grown from a polycrystalline feed 
rod employing the floating zone optical furnace. The X- 
ray characterization shows that the a- and c-axes signifi- 
cantly decrease with the Zn substitution (Fig. 1) whereas 
the b-axis shows only a minor decrease within the reso- 
lution of the measurements. The smooth change of a, 
b, and c indicates that a solid solution of MnW04 and 
ZnW04 has formed. Small pieces were cut from the crys- 
tals and oriented using Laue single crystal X-ray tech- 
niques. The size and shape of the crystals was carefully 
chosen to fit the demands of the different experiments. 
The susceptibility was measured in a 5 Tesla Magnetic 
Property Measurement System (MPMS, Quantum De- 
sign). The ferroelectric polarization of the multiferroic 
phase was measured by integrating the pyroelectric cur- 
rent upon heating the sample in zero electric field at a 
speed of 1 K/min. Before the pyroelectric measurement 
the sample was cooled in an electric bias field of 3 kV/ cm 
to align the ferroelectric domains. The heat capacity was 
measured in the Physical Property Measurement System 
(PPMS, Quantum Design). The elastic neutron scatter- 
ing measurements were performed at the HB1A and HB1 



III. RESULTS AND DISCUSSION 

The b-axis magnetic susceptibility at high tempera- 
tures (T>50 K) follows the Curie- Weiss law, as shown in 
Fig. 2. From the inverse susceptibility the effective mag- 
netic moment /i e // per formula unit (Mni_ 2; Zn 3 ;W04) 
as well as the Curie- Weiss temperature @w can be ex- 
tracted and are displayed in the inset of Fig. 2. ®w 
is negative as expected for dominantly antiferromagnetic 
coupling and it's magnitude decreases with Zn substitu- 
tion indicating the weakening of the magnetic coupling 
between the Mn spins. The effective magnetic moment 
also decreases with x since the Zn substitution contin- 
uously dilutes the system of Mn spins. The expected 
values for /i e //, assuming localized spins of the Mn 2+ 
ions with S=5/2, are shown by the dashed line in the in- 
set of Fig. 2. The experimentally determined values for 
fj, e f f are in very good agreement with the calculated val- 
ues confirming the conclusion that the Zn 2+ substitutes 
for the Mn 2+ at levels up to the maximum of this work, 
50 %. This is further confirmed by Inductively Coupled 
Plasma Mass Spectrometry, testing the elemental com- 
position at ten different spots of a single crystal with 
nominally 50 % Zn doping. The average Zn concentra- 
tion was determined as 0.49 ±0.03, close to the nominal 
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FIG. 3. (Color online) Low-temperature magnetic suscepti- 
bility Xb of Mni- I Zn I W04. Different curves are vertically 
offset for enhanced clarity. 



composition. 

The b-axis magnetic susceptibility at low temperatures 
clearly reveals the magnetic phase transitions at TV and 
Tl, as shown in Fig. 3 (different curves are vertically off- 
set for better clarity) . The transition into the AF3 phase 
is defined by the maximum of Xb(T) as indicated by the 
dashed line labeled T^- With the lock-in transition into 
the AF1 phase Xb shows a sharp increase at T^ which is 
quickly shifted to lower temperatures with increasing x 
and is completely missing for x>0.05. 5% of Zn substitu- 
tion is sufficient to suppress the collinear AF1 phase and 
to stabilize the multiferrroic ferroelectric phase as the 
ground state. The effect of Zn substitution is opposite 
to the results of Fe substitutional where the multiferroic 
AF2 phase was suppressed and the low-T AF1 phase was 
stabilized with increasing Fe content. It is also interesting 
to note that a similar stabilization of the AF2 phase was 
reported for Co substituted MnW04 powder samples^ 
although the magnetic phases of Mrii_ x Co x W04 are far 
more complex and several coexisting commensurate and 
incommensurate phase have been observed in single crys- 
tals at higher Co doping.— 

The main feature of the multiferroic phase is the ex- 
istence of a spontaneous ferroelectric polarization that 
is oriented along the b-axis. Fig. 4 shows P&(T) of 
Mni_ 2; Zn a ;W04 for x between and 0.5 as measured by 
the pyroelectric current method. In perfect agreement 
with the magnetic data, the ferroelectric phase extends 
to the lowest temperature for x values equal or above 5 %. 
The multiferroic phase becomes the ground state and it 
exists even at substitution levels as high as 50 % while the 
onset temperature, Tc, is reduced with increasing x. It 
is interesting that the 2 % substitution does decrease 
rapidly but there remains a finite polarization in the low- 
temperature state indicating the coexistence of the AF1 
and AF2 phases. Similar phase coexistence has also been 
observed in recent single crystal neutron scattering exper- 
iments on Mn^Fe^WO^^ and Mno.85Coo.15WO4.22 
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FIG. 4. (Color online) Ferroelectric polarization of 
Mni__ a; Zn a: W04 in zero magnetic field. Tat and Tl are la- 
beled for the x=0 data. 



At higher substitution levels the ferroelectric polariza- 
tion increases continuously with decreasing temperature. 
The polarization data of Fig. 4 prove that, unlike the 
Fe-substituted system Mni_ a; Fe a ;W04, the substitution 
of the non-magnetic Zn ion favors the multiferroic ferro- 
electric state which becomes the ground state for x>0.05. 

To identify the magnetic orders and their development 
with increasing Zn substitution neutron scattering ex- 
periments have been conducted. The temperature de- 
pendence of the integrated magnetic peak intensities in 
doped Mni_ a; Zn 2 ;W04 is displayed in Fig. 5. The mea- 
surements were performed near the characteristic mag- 
netic wavevector lf 2 = (-0.214, 1/2,0.457) for the AF2 
phase that is associated with the ferroelectric behavior. 
At the lowest Zn concentration x=0.02, the peak inten- 
sity first increases upon cooling below 13 K, and is sig- 
nificantly suppressed below 6 K but remains finite at 
lower temperature. Such behavior is accompanied by the 
sudden appearance of the commensurate magnetic order 
with -f 1 = (-0.25,0.5,0.5) below T L =6 K (solid circles 
in Fig. 5). For higher Zn-doping systems, only incom- 
mensurate magnetic orders exist at all temperatures, and 
the thermal evolution of the integrated intensity closely 
tracks the polarization measurement evidencing the inti- 
mate connection between the spiral order and the ferro- 
electricity. 

At the high-temperature side the onset of magnetic or- 
der (AF3 phase) is clearly visible in the sudden increase 
of the neutron intensities at T^r. The transition from 
the sinusoidal AF3 phase to the helical AF2 phase at 
Tc is more difficult to detect because both phases have 
the same magnetic modulation vector. A small change 
of slope and an anomaly in the width of the magnetic 
scattering peak is associated with the AF3— »AF2 transi- 
tion. The derived critical temperatures, Tn and Tc, are 
included in the phase diagram of Fig. 7 (triangles). 
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FIG. 5. (Color online) Neutron scattering intensities of the in- 
commensurate (open symbols, q2 = (—0.214, 1/2,0.457)) and 
commensurate (filled symbols, qi — (—0.25, 0.5, 0.5)) peaks of 
Mni- I Zn 1 W0 4 . 



Similarly, the two transitions at Tjv and Tc are also 
difficult to distinguish in the magnetization data of Fig. 3. 
The heat capacity C P (T), however, shows two well- 
defined sharp anomalies at the two phase transitions, a 
sharp rise and peak at and a second peak at Tc . We 
have measured the heat capacity for all samples as shown 
in Fig. 6. Both transitions are well resolved although the 
width of the transitions slightly increases with x. For the 
undoped M11WO4 (x=0) the sharp peak at 7.8 K indi- 
cates the first order phase transition into the AF1 phase. 
The T-x phase diagram is constructed from heat capacity, 
magnetization, and polarization measurements (Fig. 7). 
The commensurate AF1 phase is completely suppressed 
for x>0.05. The finite value of the low-temperature po- 
larization for x=0.02 reveals the coexistence of AF1 and 
AF2 below T^, as confirmed by the neutron scattering 
experiments (Fig. 5). 

The two phases (AF1 + AF2) coexisting for x<0.05 
at low T can be tuned by external magnetic fields and 
the degeneracy is completely lifted above a critical field. 
This is best demonstrated in measuring the ferroelec- 
tric polarization of Mno.gsZno.o2W04 in external mag- 
netic fields with different orientations (Fig. 8). With the 
field directed along the spin easy axis (a-c plane, Fig. 
8a) the low-T polarization (T<T^) rapidly increases and 
for H>2 T the small drop of P 6 at T L is gone. The 
magnetic field has completely suppressed the AF1 phase 
and the multiferroic AF2 phase extends entirely to zero 
temperature. However, with the field applied along the 
b-axis (Fig. 8b) the low-T polarization decreases and 
reaches zero above 3 T. The low temperature phase and 
the ground state is now the AF1 phase and the degen- 
eracy of both phases is removed. The lock-in transition 
temperature increases with from 5.7 K (zero field) to 
8.7 K at 7 T. 




FIG. 6. (Color online) Heat capacity, C p /T, of 
Mni_ a; Zn a ;W04 in zero magnetic field. Different curves are 
vertically offset by two units (zero is defined by the dashed 
line for each curve). 



IV. SUMMARY AND CONCLUSIONS 

We have shown that the multiferroic AF2 phase in 
Mni-zZn^WCU is remarkably stable with respect to the 
dilution of the magnetic exchange couplings induced by 
the substitution of nonmagnetic Zn ions. At low substi- 
tution levels and low temperatures a coexistence of two 
phases, the paraelectric AF1 phase and the ferroelectric 
AF2 phase, is shown. External magnetic fields do resolve 
the degeneracy of the two phases. Depending on the ori- 
entation of the magnetic field cither one of the coexisting 
states may become the ground state. 

The rapid suppression of the collinear AF1 phase by 
less than 5 % Zn doping and the survival of the sinu- 
soidal AF3 and helical AF2 phases up to 50 % substitu- 
tion is unique and distinguishes the multiferroic MnWCU 
from other multiferroics with a similar phase sequence, 
e.g. Ni3V2 08^ In the latter compound it was shown 
that the magnetic and multiferroic phases are strongly 
suppressed with the substitution of the magnetic Ni by 
nonmagnetic Zn, consistent with a two-dimensional spin 
system^ In MnWC>4, however, the relatively small ef- 
fect of the dilution of the system of Mn-spins through Zn 
doping on the AF3 and AF2 critical temperatures indi- 
cates a three-dimensional character of the basic magnetic 
exchange interactions. The microscopic exchange inter- 
actions can be revealed through inelastic neutron scatter- 
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FIG. 7. (Color online) Multiferroic phase diagram of 
Mni_ a; Zn a: W04. The low-temperature phase below x=0.05 
is a mixed phase, AF1 + AF2. Squares denote data from 
magnetic, heat capacity, and polarization measurements. Tri- 
angles are data from neutron scattering. 
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FIG. 8. (Color online) Magnetic field effect on the polariza- 
tion of Mno.98Zno.o2W04. The phase coexistence at low T is 
lifted in external fields. 

ing (INS) experiments probing the magnetic excitations. 
By comparing the magnon spectra with standard models 



the most relevant exchange parameters can be extracted. 
According to recent INS experiments on MnW0 4 short 
range (nearest and next-nearest neighbors) exchange in- 
teractions are not sufficient to explain the magnetic exci- 
tation spectrum but up to 11 different exchange pathways 
have to be involved to fit the data i 42 ' 43 The higher or- 
der exchange coupling constants involved in the magnetic 
correlations prove the three-dimensional character of the 
magnetic fluctuations. In three dimensions the perco- 
lation threshold for site dilution is much lower than for 
two-dimensional systems explaining the robustness of the 
magnetic orders (AF3 and AF2) in Zn-doped MnW04. 

The low-temperature AF1 phase, however, is in strong 
competition with the helical AF2 state. A simple mean- 
field model calculation has shown that the phase bound- 
ary between AF1 and AF2 as function of competing ex- 
change coupling and anisotropy constants can be very 
steep and extremely sensitive with respect to small 
perturbationsi 20 i 44 The frustration of different magnetic 
states results in the deviation from the mean field result 
for the suppression of the critical temperature as a func- 
tion of doping level in a diluted magnetic system and 
the strong suppression of the AF1 phase as observed in 
Mni_ 2; Zn 2 ,W04. 

Note added: After submission of this manuscript the 
authors became aware of a related investigation of poly- 
crystalline MnW0 4 doped with Mg and Zn up to 30 %4£ 
The conclusions derived from magnetization and dielec- 
tric constant data are consistent with the phase diagram 
of Fig. 7. 



ACKNOWLEDGMENTS 

This work is supported in part by the T.L.L. Temple 
Foundation, the J. J. and R. Moores Endowment, and the 
State of Texas through TCSUH and at LBNL through 
the US DOE, Contract No. DE-AC03-76SF00098. The 
research at Oak Ridge National Laboratory's High Flux 
Isotope Reactor was sponsored by the Scientific User Fa- 
cilities, Office of Basic Energy Sciences, U. S. Department 
of Energy. 



1 M. Fiebig, J. Phys. D: Appl. Phys. 38, R123 (2005). 

2 N. A. Spaldin and M. Fiebig, Science 309, 391 (2005). 

3 Y. Tokura, J. Mag. Mag. Mat. 310, 1145 (2007). 

4 M. Kenzelmann, A. B. Harris, S. Jonas, C. Broholm, 
J. Schefer, S. B. Kim, C. L. Zhang, S.-W. Cheong, O. P. 
Vajk, and J. W. Lynn, Phys. Rev. Lett. 95, 087206 (2005). 

5 M. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006). 

6 G. Lawes, A. B. Harris, T. Kimura, N. Rogado, R. J. Cava, 
A. Aharony, O. Entin-Wohlman, T. Yildirim, M. Kenzel- 
mann, C. Broholm, and A. P. Ramirez, Phys. Rev. Lett. 
95, 087205 (2005). 



7 K. Taniguchi, N. Abe, T. Takenobu, Y. Iwasa, and 
T. Arima, Phys. Rev. Lett. 97, 097203 (2006). 

8 O. Heyer, N. Hollmann, I. Klassen, S. Jodlauk, L. Bohaty, 
P. Becker, J. A. Mydosh, T. Lorenz, and D. Khomskii, J. 
Phys.: Condens. Matter 18, L471 (2006). 

9 H. Katsura, N. Nagaosa, and A. V. Balatsky, Phys. Rev. 
Lett. 95, 057205 (2005). 

10 I. A. Sergienko and E. Dagotto, Phys. Rev. B 73, 094434 
(2006). 

11 I. A. Sergienko, C. Sen, and E. Dagotto, Phys. Rev. Lett. 
97, 227204 (2006). 



6 



12 M. Mochizuki, N. Furukawa, and N. Nagaosa, Phys. Rev. 
Lett. 105, 037205 (2010). 

13 D. Higashiyama, S. Miyasaka, N. Kida, T. Arima, and 
Y. Tokura, Phys. Rev. B 70, 174405 (2004). 

14 N. Hur, S. Park, P. A. Sharma, S. Guha, and S.-W. 
Cheong, Phys. Rev. Lett. 93, 107207 (2004). 

15 S. Scki, Y. Yamasaki, M. Soda, M. Matsuura, K. Hirota, 
and Y. Tokura, Phys. Rev. Lett. 100, 127201 (2008). 

16 C. R. dela Cruz, B. Lorenz, Y. Y. Sun, Y. Wang, S. Park, 
S.-W. Cheong, M. M. Gospodinov, and C. W. Chu, Phys. 
Rev. B 76, 174106 (2007). 

17 C. R. dela Cruz, B. Lorenz, and C. W. Chu, Physica B 
403, 1331 (2008). 

18 R. P. Chaudhury, F. Yen, C. R. dela Cruz, B. Lorenz, Y. Q. 
Wang, Y. Y. Sun, and C. W. Chu, Phys. Rev. B 75, 012407 

(2007) . 

19 R. P. Chaudhury, C. R. dela Cruz, B. Lorenz, Y. Sun, 
C. W. Chu, S. Park, and Sang-W. Cheong, Phys. Rev. B 
77, 220104(R) (2008). 

20 R. P. Chaudhury, B. Lorenz, Y. Q. Wang, Y. Y. Sun, and 
C. W. Chu, Phys. Rev. B 77, 104406 (2008). 

21 R. P. Chaudhury, B. Lorenz, Y.-Q. Wang, Y. Y. Sun, and 
C. W. Chu, New J. Phys. 11, 033036 (2009). 

22 S. Scki, Y. Yamasaki, Y. Shiomi, S. Iguchi, Y. Onose, and 
Y. Tokura, Phys. Rev. B 75, 100403(R) (2007). 

23 S. Kanetsuki, S. Mitsuda, T. Nakajima, D. Anazawa, H. A. 
Katori, and K. Prokes, J. Phys.: Condons. Matter 19, 
145244 (2007). 

24 G. Lautenschlager, H. Weitzel, T. Vogt, R. Hock, A. Bohm, 
M. Bonnet, and H. Fuess, Phys. Rev. B 48, 6087 (1993). 

25 H. Sagayama, K. Taniguchi, N. Abe, T. H. Arima, 
M. Soda, M. Matsuura, and K. Hirota, Phys. Rev. B 77, 
220407(R) (2008). 

26 H. Ehrenberg, H. Weitzel, C. Heidy, H. Fuess, G. Wltschek, 
T. Kroener, J. van Tol, and M. Bonner, J. Phys. Condens. 
Matter 9, 3189 (1997). 

27 A. H. Arkenbout, T. T. M. Palstra, T. Siegrist, and 
T. Kimura, Phys. Rev. B 74, 184431 (2006). 

28 K. Taniguchi, N. Abe, H. Sagayama, S. Ohtani, 
T. Takenobu, Y. Iwasa, and T. Arima, Phys. Rev. B 77, 
064408 (2008). 

29 R. P. Chaudhury, F. Yen, C. R. dela Cruz, B. Lorenz, Y. Q. 
Wang, Y. Y. Sun, and C. W. Chu, Physica B 403, 1428 

(2008) . 



30 D. Meier, M. Maringer, T. Lottermoser, P. Becker, L. Bo- 
haty, and M. Fiebig, Phys. Rev. Lett. 102, 107202 (2009). 

31 T. Finger, D. Senff, K. Schmalzl, W. Schmidt, L. P. Reg- 
nault, P. Becker, L. Bohaty, and M. Braden, Phys. Rev. B 
81, 054430 (2010). 

32 H. A. Obermayer, H. Dachs, and H. Schrocke, Solid State 
Commun. 12, 779 (1973). 

33 C. Klein and R. Geller, J. de Physique 35, C6:589 (1974). 

34 E. Garcia-Matres, N. Stiifier, M. Hofmann, and 
M. Reehuis, Eur. Phys. J. B 32, 35 (2003). 

35 Y. Ding, N. Stiifier, M. Reehuis, M. Hofmann, G. Ehlers, 
D. Giinther, M. Meifiner, S. Welzel, M. Wilhelm, 
M. Steiner, and F. Kubanek, Physica B 276-278, 596 
(2000). 

36 Y.-S. Song, J.-H. Chung, J. M. S. Park, and Y.-N. Choi, 
Phys. Rev. B 79, 224415 (2009). 

37 F. Ye, Y. Ren, J. A. Fernandez-Baca, H. A. Mook, J. W. 
Lynn, R. P. Chaudhury, Y.-Q. Wang, B. Lorenz, and C. W. 
Chu, Phys. Rev. B 78, 193101 (2008). 

38 G. Lawes, M. Kenzelmann, N. Rogado, K. H. Kim, G. A. 
Jorge, R. J. Cava, A. Aharony, O. Entin-Wohlman, A. B. 
Harris, T. Yildirim, Q. Z. Huang, S. Park, C. Broholm, 
and A. P. Ramirez, Phys. Rev. Lett. 93, 247201 (2004). 

39 R. P. Chaudhury, F. Ye, J. A. Fernandez-Baca, Y. Q. 
Wang, Y. Y. Sun, B. Lorenz, H. A. Mook, and C. W. 
Chu, Phys. Rev. B 82, 184422 (2010). 

40 R. P. Chaudhury, B. Lorenz, Y. Q. Wang, Y. Y. Sun, C. W. 
Chu, F. Yc, J. Fernandez-Baca, H. Mook, and J. Lynn, J. 
Appl. Phys. 105, 07D913 (2009). 

41 P. Kharel, A. Kumarasiri, A. Dixit, N. Rogado, R. J. Cava, 
and G. Lawes, Phil. Mag. 89, 1923 (2009). 

42 H. Ehrenberg, H. Weitzel, H. Fuess, and B. Hennion, J. 
Phys.: Condens. Matter 11, 2649 (1999). 

43 F. Yc, R. S. Fishman, J. A. Fernandez-Baca, A. A. 
Podlesnyak, G. Ehlers, H. A. Mook, Y.-Q. Wang, B. 
Lorenz, C. W. Chu, submitted for publication, Dec. 1, 
2010. 

44 M. Kenzelmann, A. B. Harris, A. Aharony, O. Entin- 
Wohlman, T. Yildirim, Q. Huang, S. Park, G. Lawes, 
C. Broholm, N. Rogado, R. J. Cava, K. H. Kim, G. Jorge, 
and A. P. Ramirez, Phys. Rev. B 74, 014429 (2006). 

45 L. Meddar, M. Josse, P. Deniard, C. La, G. Andre, 
F. Damay, V. Petricek, S. Jobic, M.-H. Whangbo, 
M. Maglione, and C. Payen, Chem. Mat. 21, 5203 (2009). 



